Lattice dynamics and structural stability of ordered Fe 3 Ni, Fe 3 Pd and Fe 3 Pt alloys 

using density functional theory 
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We investigate the binding surface along the Bain path and phonon dispersion relations for the 
cubic phase of the ferromagnetic binary alloys Fe^X (X — Ni, Pd, Pt) for LI2 and DO22 ordered 
phases from first principles by means of density functional theory. The phonon dispersion relations 
exhibit a softening of the transverse acoustic mode at the M-point in the Ll2-phase in accordance 
with experiments for ordered FesPt. This instability can be associated with a rotational movement 
of the Fe-atoms around the Ni-group element in the neighboring layers and is accompanied by an 
extensive reconstruction of the Fermi surface. In addition, we find an incomplete softening in [111] 
direction which is strongest for FeaNi. We conclude that besides the valence electron density also 
the specific Fe-content and the masses of the alloying partners should be considered as parameters 
for the design of Fe-based functional magentic materials. 
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I. INTRODUCTION 



Ferrous alloys are well known for a wide variety of 
anomalous structural and magnetic properties. Around 
1185 K, bulk iron undergoes a structural transition from 
a face centered cubic (fee) high-tcmpcraturc phase to a 
low-tempcrature body centered cubic (bec) phase, which 
is ferromagentic with a Curie temperature of 1043 K. 
The structural transition temperature is systematically 
lowered by alloying elements from the right side of Fe 
in the periodic table, which effectively lowers the va- 
lence electron concentration e/a. The transition is dis- 
placive and diffusionless between a high-tcmpcraturc 
high-symmetry austenitic phase and a low-tempcraturc 
martensitic phase, in which the symmetry is usually low- 
ered, e.g., by tetragonal distortions. In ferromagnetic 
Fe-alloys with e/a between 8.7 and 8.5, the marten- 
sitic transition disappears and fee austenite becomes 
the stable ground state structure. Iron alloys with Ni- 
group elements at this valence electron concentration 
are at the center of scientific interest for a long time, 
as they show, in addition to the martensitic instabil- 
ity, a multitude of magneto-structural anomalies that 
are of technological importance. The probably most 
prominent is the so-called Invar effect, describing the 
over-compensation of thermal expansion in a wide tem- 
perature range, which has been intensively studied for 
more than one century* 1 - - — The best known representa- 
tive is Fe65Ni35, but similar behavior is also observed 
in Fe7oPd3o and Fc75Pt25. Although most explanations 
proposed so far rely on magneto-volume coupling, the 
vicinity to the martensitic transformation has been well 
noticed and discussed. The relevant literature on this 
topic is too numerous to be summarized at this point - 
for an introduction and further references, see, e. g., Refs. 

m 



During the past three decades, the so-called magnetic 
shape memory (MSM) effect, which is again present in 
the alloys with all three Ni-group elements, was brought 
to the attention of the scientific community. The MSM 
behavior allows very large magnctic-ficld-induced strains 
of up to three percent to be achieved by externally apply- 
ing moderate magnetic fields in the (sub-)Tesla ranged— 
The effect is not as pronounced as for the Ni-based 
Heusler alloys like Ni 2 MnGa^-^ but its origin is be- 
lieved to be related. One common key ingredient is the 
appearance of intermediate or modulated martensites - 
5-fold and 10-fold modulated pseudo-tetragonal or ortho- 
hombic in the case of Ni-Mn-Ga^— and (slightly dis- 
torted) face centered tetragonal (fct) in the case of Fc-Pt 
and Fe-Pd - between the high symmetry austenite and 
the low-temperature martensite, which is body centered 
tetragonal in the case of the ferrous alloysJ^— 

For the above mentioned alloys, the relevant electron 
concentration range is achieved at compositions which 
may allow the formation of ordered stoichiometric com- 
pounds. However, FeaPt is the only alloy where the 
ordered phase is rcproducibly realized in experiment. 
Nevertheless, measurements are predominately made for 
slightly off-stoichiometric ordered Fe72Pt28 which does 
not transform martensitically, even if order is not fully 
complete. Therefore, most data are available only for 
disordered or partially ordered alloys. Despite this fact, 
we will restrict our study to stoichiometric, ordered com- 
pounds. These are much easier to handle from the techni- 
cal point of view and allow deeper insight into electronic 
and vibrational properties with respect to their physical 
origin, because their simulation cells are sufficiently small 
while side-effects as statistical broadening can be largely 
avoided. 

Another important feature of these alloys is the marked 
softening of the transversal acoustic phonons in [110] di- 
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FIG. 1: Primitive cells of fee LI2 (a) and DO22 (b) ordered 
alloys of A3B stoichiometry as used in our calculations. The 
lower two images show both structures represented in the 
same supercell which consists of two primitive cells stacked 
along the c-axis in the case of L2i (c). Exchanging the A 
and B components (marked by dark blue and bright orange 
spheres) in one of the a-6-planes yields the DO22 ordering (d). 
Note, that the symmetry of the DO22 structure is tetragonal. 
After reduction of the c/a-ratio by a factor of a/1/2, the cell 
transforms to the bec-type DO3 structure with cubic symme- 
try. 



rection, which has subsequently been related to any of 
the three above mentioned anomalies in the pasti^ir— 
Our central aim is thus to investigate the influence of 
electron-phonon coupling in FeaNi, FeaPd and FeaPt as 
one of the important mechanisms for structural trans- 
formations of Fe-based binary alloys. This complements 
the experimental information by providing a link between 
electronic properties and structural distortions, which 
will help to clarify the microscopic origin of these insta- 
bilities. We accomplished this aim by comparing features 
of the phonon dispersion with the electronic structure of 
these alloys. A similar approach has been used in the 
past to explore the nature of electron-phonon coupling 
in conventional and magnetic shape memory alloys^r— 
During our preoccupation with this topic, we noticed 
that investigations providing a structural and dynamical 
view of the field on the basis of first principles calcu- 
lations are sparse and performed with varying method- 
ology, hindering a comparison of all three isoelectronic 
alloys on equal footing. We therefore decided to provide 
first a systematic survey of structural properties and lat- 
tice dynamics of ordered Fe-rich alloys with elements of 
the Ni-group, thereby filling the above mentioned gap 
in existing literature. The experimentally observed or- 
der in Fc 3 Pt is of Ll 2 type, see Fig. [TJ In our calcu- 
lations, we will alternatively consider DO22 order. For 
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FIG. 2: The energy profile along the Bain path for LI2 
Fe 3 Pd obtained with VASP and Wien2k for PBE and PW9I 
exchange-correlation potentials. The calculations were car- 
ried out at a fixed atomic volume of 13.08 A 3 , corresponding 
to an fee lattice constant of a = 3.74 A. For better compari- 
son, all results are given relative to the energy of the bee state 
(c/a = a/1/2). The labeling fee and bee refers to coordina- 
tion of the atoms regardless of the species. With the VASP 
code, pseudopotentials with different numbers of semi-core 
electrons were tested; the reference calculations (Wien2k) ac- 
count for Fe 3s and Pd 4s as semi-core electrons. Most curves 
coincide along the full Bain path within 2meV/atom. 

lattices with bec coordination DO22 becomes equivalent 
with DO3 order, which is related to the L2i structure of 
ternary Heusler alloys. 

Fe-rich Fe-Pd and Fe-Pt alloys close to the stoichio- 
metric composition are reported to be collinear ferromag- 
nets at low temperatures; we will therefore solely concen- 
trate on the collinear ferromagnetic case. Excited non- 
collinear, ferri- and antifcrromagnetic spin structures, 
which have been discussed with respect to the Invar 
effect r^r— would certainly be beneficial for a thorough 
understanding of the interdependence of lattice dynamics 
and finite temperature anomalies as MSM and Invar ef- 
fect. These, however, arc beyond the scope of the current 
work and are thus left open for future investigation. 

II. NUMERICAL DETAILS 

Our investigation of the structural and dynamic prop- 
erties has been performed within the framework of den- 
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sity functional theory (DFT)/ 3 ^ The majority of our re- 
sults has been obtained using the Vienna ab initio Sim- 
ulation package (VASP) . 39 ' 40 This code concentrates - 
as pseudopotential methods generally do - for efficiency 
reasons on the valence electrons for the description of 
the electronic structure. The VASP code provides, nev- 
ertheless, an excellent compromise between speed and ac- 
curacy as the interaction with core electrons is taken into 
account within the projector augmented wave approach 
(PAW).— This yields results close to what is usually ex- 
pected from other all electron methods. As we are deal- 
ing with Fe-rich compositions, the use of the generalized 
gradient approximation (GGA) is mandatory to obtain a 
correct description of the structural ground state prop- 
erties. For structural energies within small (four atom) 
cells corresponding to Fig. [T] (a) and (b) we use a fc- 
mesh containing at least 16384 fc-points x atoms in the 
full Brillouin zone. For the calculation of the density of 
states (DOS) the fc-point density was increased to values 
of 131072 fc-points x atoms and 340736 fc-points x atoms 
for obtaining the Fermi surface in the 4 atom cell. The 
Brillouin zone integration was carried out using the tetra- 
hedron method with Blochl corrections^ Electronic self- 
consistency was assumed below a threshold difference of 
1 iieV between two iteration. The VASP code is supplied 
with an extensive, well tested potential library which con- 
tains more than one pre-generated pseudopotential per 
element. These differ by number of electrons explicitly 
treated as valence or semi-core states or are designed for 
different exchange correlation functionals. In our case, 
the GGA formulations of Perdew and Wang (PW91) in 
connection with the spin-interpolation formula of Vosko, 
Wilk and Nusair as well as Perdew Burke and Ernzer- 
hof (PBE) have been used.— ~— The required cutoff of 
the plane wave basis depends on the exchange corre- 
lation potential and was chosen as E cut = 335 eV (or 
larger) for PW91 and PBE together with the minimal 
basis explicitly describing the 3c? 7 4s 1 for Fe and 4d 9 5s 1 
for Pd (and correspondingly for Ni and Pt). In calcu- 
lations cxplicitely considering the semi-core p-electrons, 
E cut was set to 367 cV for PBE and 438 cV for the PW91 
pseudopotentials. A scalar relativistic formulation of the 
Hamiltonian was used throughout 4£ 



A. Comparison between VASP and Wien2k 

Although the above mentioned technical differences 
with respect to the choice of the exchange correlation 
functional and the basis size usually do not lead to qual- 
itative changes of the results, we must be particularly 
careful in our case, because energy landscapes in the 
vicinity of a structural transformation can be essentially 
flat. This requires a higher accuracy, i. e., a better energy 
resolution. Therefore, we decided to perform in a first 
step a thorough comparison of the available potentials. 
As a benchmark, we evaluated the energy as a function 
of the tetragonal distortion c/a along the Bain path 4 ^ for 



LI2 ordered FcaPd, describing a structural deformation 
from a face centered cubic (fee) lattice at c/a = 1 to 
a body centered cubic (bec) lattice at c/a = \Jl/2, as 
shown in Fig. [2j The results were compared to high pre- 
cision calculations using the full potential linearized aug- 
mented plane wave method (FLAPW) as implemented 
within Wien2k code^ 9 - which is widely counted among 
the most accurate DFT codes for solid state problems. 
Here, the calculations included 3s and 3p of iron and 4s 
and Ap of palladium as semi-core. The energetic cut- 
off was chosen as Rmt x K max = 9.0. Muffin-tin radii 
of 2.24 a.u. were used for all atoms and angular mo- 
menta were taken into account up to l max = 10. The 
fc-mesh comprised again 32000 fc-points x atoms in the ir- 
reducible Brillouin zone. 

The obtained energy variation along the Bain path is 
with m 8meV/atom very small. However, except for the 
PBE potentials without semi-core p electrons, all curves 
coincide within an energy interval of 2. . .3 meV/atom 
(23 ... 35 K on a temperature scale) , which can be taken 
as a measure for the methodological resolution in calcu- 
lating structural energy differences along the Bain path. 
While the PBE potentials seem to systematically over- 
estimate the energy variation (slightly), the agreement 
between both PW91 potentials (large and small basis) 
and the full potential curves is very good. This finding 
agrees well with a recent comparison of experimental and 
theoretical structural data for equiatomic FePt alloys^ 
Our benchmark demonstrates that the PW91 pseudopo- 
tentials with the small basis (Fe 3p and Pd Ap> as core 
electrons) is sufficient, which is especially helpful for ex- 
pensive calculations of large systems. 



B. Calculation of dynamical properties 

Two standard techniques are currently applied for the 
investigation of lattice dynamics of crystals from first 
principles: The linear response method and the so-called 
direct method^"— 

In the linear response method the dynamical matrix 
is obtained from the modification of the electron den- 
sity, via the inverse of dielectric matrix describing the re- 
sponse of the valence electron density to a periodic lattice 
perturbation. The dielectric matrix is then calculated 
from the eigenfunctions and energy levels of the unper- 
turbed system^ Only linear effects, such as harmonic 
phonons, are accessible to this technique. The method 
has been applied with success to many alloys related to 
our present stud y 31 ' 55 

On the other hand, the direct-method is a frozen- 
phonon type of calculations based on a supercell calcula- 
tion, which allows explicit account of any distortion of the 
atomic positions. The phonon frequencies are calculated 
from Hellmann-Feynman forces generated by the small 
atomic displacements, one at a time. Hence using the in- 
formation of the crystal symmetry space group, the force 
constants are derived, the dynamical matrix is built and 
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diagonalized, and its eigenvalues arranged into phonon 
dispersion relations. In this way, phonon frequencies at 
selected high-symmetry points of the Brillouin zone can 
be calculated^ The direct approach in conjunction with 
the aft-initio method has previously been used intensively 
by the authors to investigate phonon dispersion relations 
in magnetic shape memory Heusler compounds! 32 ' 57 " — 

The results presented within this work are relying on 
the direct approach together with the VASP code for 
the calculation of the respective forces; the displacements 
being necessary to describe the phonon dispersions were 
generated by the PHON code written by Dario Alfe i 60 ' 61 
This code was also used later on to generate the dynami- 
cal matrix and the resulting phonon dispersion relations. 
We used a supercell of 5x5x5 primitive cells, containing 
500 atoms in total. For sufficient accuracy of the forces 
a fc-mesh of 4x4x4 points in the full Brillouin zone was 
employed in connection with the Methfessel-Paxton fi- 
nite temperature integration scheme (smearing param- 
eter er = 0.1eV). First-principles calculations of this size 
are computationally very demanding and can, so far, only 
be performed on world leading supercomputer installa- 
tions. The forces were calculated for three independent 
displacements of 0.02 A in size each. From these calcula- 
tions the phonon dispersions along the lines connecting 
the main symmetry points and the vibrational density of 
states (VDOS) were calculated. For the latter, a mesh of 
61x61x61 points in reciprocal q-space was used and ad- 
ditional Gaussian broadening with a smearing parameter 
of <r = 0.1 THz was applied. 

We also carried out comparative calculations for DO3 
FeaNi and LI2 Fe 3 Pt using the linear response method 
as implemented in the PWSCF packaged To describe 
the interaction between ionic cores and valence elec- 
trons, we used for Fe3Ni ultra-soft pseudo-potentials^ 
generated using the exchange correlation functional of 
Perdew, Burke and Ernzerhof (PBE). The pseudopoten- 
tials for Fe 3 Pt were generated for the exchange correla- 
tion of Pcrdcw-Zunger within local density approxima- 
tion (LDA)M- The technical parameters used here, were 
a kinetic energy cutoff of 50 Ry, an energy cutoff for 
the augmentation charges of 600 Ry and 12x12x12 k- 
points in the full Brillouin zone (Monkhorst-Pack) . For 
the phonon dispersions, a mesh of 4x4x4 q-points in the 
reciprocal space was used for the fee basis, yielding 8 
sets with finite weight while for the simple cubic basis 
a 2x2x2 q-points mesh was used yielding 4 sets of spe- 
cial q- vectors with finite weight. Although the technical 
parameters are less tight, the obtained dispersions agree 
well on a qualitative level with the ones obtained within 
the direct approach (except for discrepancies which can 
be partially related to the use of LDA in the latter case 
as the typical underestimation of the lattice constant as 
well as an extension of the instability of the TAi branch 
to the T-point, see also Ref. [65h and are therefore not 
explicitly presented in this manuscript. 



III. COMPARISON OF STRUCTURAL 
PROPERTIES 

There are a large number of reports of first principle 
investigations concerning Fc-rich alloys with elements of 
the platinum group available in literature. These cover 
magnetic and magneto-elastic properties* 6 ^— electronic 
structure ; 69 ' 73 " — order and disorder i 7 - 6 -"— Only few re- 
ports, however, are present describing the energetics of 
tetragonal distortions in ordered alloys. These exist, e. g., 
for Fe 3 Pt using the LMTO-ASA (linear muffin-tin or- 
bital method within the atomic sphere approximation) 
method and for FeaNi using different approaches! 69 ' 79 " — 
Again, particularly the FLAPW calculations of Ref. [sT| 
agree very well with our results. A systematic ab-initio 
comparison of all three alloys with respect to structural 
properties along the Bain path, describing the transfor- 
mation from fec-type to bec-type structure, however, is 
still missing and therefore concisely presented in the sub- 
sequent section. 



A. Binding surfaces along the Bain path 

The total energy as a function of atomic volume and 
tetragonal distortion for all three isoelectronic alloys in 
LI2 order is shown in Fig. [3J As is easily seen from the 
contour plots, the optimum c/a ratio varies from close 
to bec (c/a = 0.74) for Fe 3 Ni over bet (c/a = 0.78) 
for FesPd to fee (c/a = 1) for Fc 3 Pt. This trend can- 
not simply be interpreted in terms of the atomic vol- 
ume as, especially for Fe 3 Pd and Fe 3 Pt, the equilib- 
rium volumes nearly coincide. This is in good agreement 
with experimental measurements (see the collection of 
Okamoto, Ref. HU, for an extensive overview). Especially 
the relation between lattice parameter and composition 
of Fe 3 Pd exhibits a strong positive deviation from Veg- 
ards law. The overall variation of the energy landscape 
along the Bain path is exceptionally small for all three 
alloys. It is with 22mcV/atom largest for Fc 3 Ni and 
smallest for Fc 3 Pd with only 4meV/atom (minimum to 
saddle point). The shift of the minima agrees well with 
the fact that in the phase diagram of ordered Fe 3 Pt the 
stability range of the body centered structures at T = K 
ends at higher Fe-content compared to the other alloys as 
well as with the trend in the martensitic transition tem- 
peratures, which are reduced on variation of the second 
element from Ni to Pt. The maximum change in equilib- 
rium volume associated with a full tetragonal transforma- 
tion along the Bain path is small, less than one percent 
for all three alloys. There is also a strong experimental 
indication for an unusually flat energy surface. Recently, 
it has been shown that thin disordered Fe7oPd 3 o films 
can be grown epitaxially on different substrates inducing 
different lattice constants a in the film plane.— The cor- 
responding perpendicular lattice constant, c, however, is 
free to adjust, thereby realizing c/a ratios covering most 
of the Bain path. Such a strained epitaxial growth can 
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FIG. 3: (color online) Contour plots of the binding sur- 
faces (energy as a function of atomic volume v and tetrago- 
nal distortion c/a) of ferromagnetic LI2 ordered FeaNi (top), 
Fe 3 Pd (center) and Fe 3 Pt (bottom). The energy distance 
between two contour lines is 2meV/atom. The global min- 
ima are located at c/a = 0.74 and v = 11.633 A 3 /atom for 
Fe 3 Ni, c/a = 0.78 and v = 13.054 A 3 /atom for Fe 3 Pd, and 
c/a = 1.00 and v = 13.082 A 3 /atom for Fe 3 Pt. 



only be expected if the energy associated with the lattice 
strain is sufficiently small along the Bain path. 

The magnetic moment is largest on the bcc side, how- 
ever, there is only little variation of the spin moments 
along the Bain path, which is of the order of 0.1 /ZB/atom 
(cf., Fig. [4j. The slight variation of the magnitude be- 
tween the alloys might well be attributed to the different 
equilibrium lattice constants. 



B. LI2 versus DO22 order 

For fee alloys of stoichiometry A3B, the DO22 structure 
may be considered as another, possibly competing, real- 
ization of order. As can be seen from Fig. [T] (d), the DO22 
structure emerges from the LI2 structure by shifting one 
of the mixed planes half way along the diagonal with- 
out shearing the other planes. At the c/a-ratio of the 
bcc lattice, the DO22 structure turns into the highly sym- 
metric DO3 structure, which is equivalent to the Heusler- 
type L2i structure for a binary composition. For FeaNi 
and FeaPd, where order is not substantiated in exper- 
iment, the DO22 becomes favored at c/a ratios slightly 
above one, as demonstrated in the lower panels of Fig. 
31 For FesPt the DO22 structure becomes favored at the 
bcc end. Only for this alloy, the energy of DO22 becomes 
competitive with L2i, taking into account the complete 
Bain path. The formation of a DO3 phase in ordered 
near-stoichiometric FcsPt might further be hindered by 
the slow kinetics of the necessary diffusive transformation 
at comparatively low temperatures. However, at the fee 
side, the small energetic distance between LI2 and DO22 
order indicates a low penalty and thus a comparatively 
high probability for anti-phase boundaries in the ordered 
phases. 

The magnetic response of both order types to a tetrag- 
onal distortion is identical. An important difference be- 
tween DO22 and L2i order is, however, that along the 
Bain path a considerable energy barrier exists for DO22 
but not for L2i, which would slow down the transforma- 
tion kinetics, if DO22 could be stabilized. 



IV. LATTICE DYNAMICS AND SOFT 
PHONONS 

Important information about the stability of structures 
and potential relaxation mechanisms can be obtained 
from the phonon dispersion relations. Accordingly, ex- 
tensive discussions have been taking place concerning 
martensitic precursor effects in the lattice dynamics of 
the austenitic phase in martensitically transforming al- 
loys, which might yield, e.g., additional information on 
the transformation mechanisms. For the systems under 
consideration, experimentally determined phonon disper- 
sion relations obtained by inelastic neutron scattering are 
available for disordered Fe7oNi 30 and Fe65Ni 35 1 29 ' 84 dis- 
ordered Fe72Pd28r^ and disordered as well as ordered 
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FIG. 4: (color online) Comparison of total energy and magnetic moment as a function of tetragonal distortion c/a between 
Ll2-ordered and D022-ordered ferromagnetic FegNi (left), Fe3Pd (center) and FesPt (right). The calculations were performed 
at a constant atomic volume close to the respective equilibrium volume of the LI2 structure. 



FeaPti 21 ' 23 " — Theoretical calculations of the full phonon 
dispersion, on the other hand, have been reported so far 
mainly for ordered Fe 3 Ni^iSL, disordered Fe 72 Pd28r^£ 
and Fe-Pd alloys in the Pd-rich composition ranged - — 
Several of the above mentioned investigations rely on em- 
pirical or semi-empirical descriptions of the interatomic 
forces. Astonishingly, similar calculations for the Fe-Pt 
case appear to be missing. The inability of the exper- 
iment to provide a complete comparison of the ordered 
phases (which are easier to understand) together with the 
incomplete record of theoretical studies motivates our at- 
tempt to provide a comparative overview of the phonon 
dispersion of all three ordered isoelectronic alloys from 
first principles. 



A. Phonon dispersion relations 

The central feature of the phonon dispersions is the 
complete softening of the transversal acoustic branch 
TAi around the M-point. This is common for all three 
alloys, but most pronounced in FeaPt. The occurrence 
of imaginary frequencies reflects the the instability of the 
lattice which may gain energy by static atomic displace- 
ments according to the corresponding phonon mode. A 
strong anomalous softening of the TAi branch is a widely 
discussed observation from neutron diffraction experi- 
ments on ordered FcsPt i 21 ! 24 ' 27 ' 28 With decreasing tem- 



perature starting from the paramagnetic phase, the TAi 
mode becomes increasingly softer at the M-point, finally 
reaching phonon energies as low as 3meV for T = 12 K. 
The microscopic nature of the corresponding displace- 
ments will be discussed in the following section. An in- 
complete softening of the TAi branch with anomalous 
temperature dependence is also found for the disordered 
alloys. However, in these cases, it is far less pronounced 
and the wave vector is shifted towards the r-pointi 22 ' 27 i 29 

For FesPd one of the transversal acoustic modes in Fig. 
[5] starts with a negative slope from the T-point, which 
indicates that the system might immediately undergo 
a long-wavelength distortion without activation barrier. 
This is not the case for FeaNi, although here, the under- 
lying structure is not the groundstate, either. In both 
cases, the cubic Ll 2 state describes only a saddle point. 
Thus the initial slope of the acoustic branches in the low 
energy and low wavevector regime depends on small de- 
tails of the binding surfaces and might well be affected 
by the specific choice of the method and technical setup, 
as demonstrated in Fig. [2] The softening along the [111] 
direction is strongest for FesNi. Both TA branches are 
degenerate, reaching frequencies as low as 1 meV around 
£ = 0.14. The corresponding eigenvectors are close to 
(1,-1,0) and (1,1,-2) and practically uniform for all 
four atoms in the primitive cells. This may correspond 
to the primary [112] shear on (111) planes, common for 
martensitic transformation mechanisms of Fe-based al- 
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FIG. 5: (color online) Calculated phonon dispersions (right panels) of LI2 ordered FesNi (top), FesPd (center) and FesPt 
(bottom) along the main symmetry directions and the corresponding vibrational density of states (left panels). Imaginary 
frequencies are omitted from the density of states. The polarization of the phonon branches according to the atomic contributions 
to the respective eigenvectors is shown in the intensity (color) coding. 
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FIG. 6: (color online) Calculated phonon dispersions (right panel) of DO3 ordered FeaNi along the main symmetry directions 
and the corresponding vibrational density of states (left panel). Again, the polarization of the phonon branches according to 
the atomic contributions to the respective eigenvectors is shown in the intensity (color) coding. 



loys (e.g., see Refs. [litED)- 

The phonon dispersion of FC3N1 as shown in the upper 
panel of Fig. [5] agrees well with previous ab initio re- 
sults for the same material obtained by the same method 
(but less restrictive accuracy)^ Differences are mainly 
present at the M-point, where the complete softening of 
the TAi mode is less pronounced in the present study and 
slightly shifted away from the Brillouin zone boundary. 
In addition, the softening in [111] direction is more pro- 
nounced in our study. Previous semi-empirical calcula- 
tions based on a tight-binding scheme on the other hand 
did only show comparable softening if the strength of 
the electron phonon coupling was artificially enhanced^ 
This underlines the necessity to fully take into account 
the coupling between electronic and lattice degrees of 
freedom within a first-principles approach. Still, the con- 
vergence with respect to the subdivision in reciprocal 
space for the electronic structure as well as for the lattice 
dynamics remains an important issue. 

Another quantity which allows for a detailed compar- 
ison with experiment is the element resolved vibrational 
density of states (VDOS) which is shown in the diagrams 
on the left of the dispersions. For alloys with a large mass 
difference, we expect a clear separation between the ele- 
mental contributions. This is indeed the case for FcaPt, 
where we find a pseudo-gap at phonon energies around 
17meV. Above this energy, there is only little contribu- 
tion from Pt vibrations, while we find a hybridization 
of Pt and Fe contributions in the energy range below. 
In FesPd, the Pd states move to higher energies, filling 
up the pseudo-gap and there is increased hybridization 
of the partial contributions at the optical bands in the 



uppermost part. This trend consequently continues for 
FeaNi, moving the central peak of the partial VDOS of 
Ni up to an energy of 22.5 meV. In turn, the distinct 
low energy peak appearing around lOmeV for FesPt and 
Fc3Pd has disappeared. In all three cases, there exists a 
significant contribution of Fe-states in the low energy re- 
gion (< 8 meV). The main peaks of the calculated VDOS 
agree well with experimental measurements of disordered 
Fe65Ni35 of Delaire and coworkers^ as well as on ordered 
FesPt by Wiele et al£^&^ and even recent measurements 
of the Fe-partial contributions of ordered FesPt nanopar- 
ticles by Roldan-Cuenya et al~± using the Mofibauer ap- 
proach. The pseudo-gap, however, is not showing up 
in experimental VDOS as these generally experience a 
stronger broadening, which may be attributed to finite 
temperature of the measurements, resolution of the mea- 
surement devices and - not at last - the incomplete or- 
dering usually present in FcsPt samples. 

The disappearance of the sharp energy peak in the 
VDOS around lOmeV points out a qualitative differ- 
ence between FeaNi and the other two alloys. In the 
latter case, these states are connected with vibrational 
modes at the M and the R-point which arc dominated 
by motions of the Pt or, respectively, Pd atoms as can 
be seen from the atomic polarization of the phonon dis- 
persions shown in the intensity (color) coding in Fig. [5] 
In FeaNi, however, the low lying modes in the vicinity 
of the M-point mainly involve Fe-atoms, which makes 
the Fe-subsystem particularly soft and susceptible to ele- 
ment specific distortions, e.g., due to magnetoelastic cou- 
pling. In contrast, this feature completely disappears in 
the body centered cubic phase of FesNi, which we rep- 
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FIG. 7: (color online) Schematic representation of the two 
M-point soft modes in the phonon dispersions of fee FesNi, 
FeaPd and FeaPt shown in Fig. [5] The picture shows a pro- 
jection of the repeated LI2 unit cell along the [0 1] direction. 
The brighter (red) arrows correspond to the M2 mode, while 
the black arrows show vibrations of the M4 mode. The main 
feature here is that both vibrations are identical, but orthog- 
onal to each other. Both vibrations match with a translation 
along one half of the the face diagonal ([110] direction) as 
shown by the thick (blue) arrow, whereby the two vibrations 
remain the same but shifted in phase by n. 
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FIG. 8: (color online) Total energy as a function of the dis- 
placement amplitude according to the two modes (M2, dia- 
monds, and M4, circles) which may be made responsible for 
the for the (complete) softening of the acoustical phonons at 
the M-point (cf. Fig. [7]). Results are shown for FesNi, FeaPd 
and FeaPt. The lines are polynomials fitted to the calculated 
data points. For all three alloys, displacements according to 
the M2-mode lead to a decrease in energy. 



resent in our calculations for simplicity by the DO3 or- 
dered structure (cf. Fig. [6]), as this also possesses cu- 
bic symmetry. Here, again, we find a rather clear sep- 
aration between the Fe and the Ni-contributions which 
now predominately occupy the low lying branches at the 
Brillouin zone boundary. In this sense, the situation 
bears similarity to the inversion of the Ni and Ga modes 
which was reported for the magnetic shape memory alloy 
Ni 2 MnGa..^2k2! 



B. Identification of the soft mode 

As laid out in the previous paragraph, the two lowest 
frequencies in the dispersion of LI2 FeaNi at the M-point 
involve only motions of the Fe species. Since the as- 
sociated wave-vector corresponds to the Brillouin zone 
boundary, the amplitudes are opposite at either end of 
the primitive cell. Therefore the term antiferrodistortive 
transformation has been has been coined for the conden- 
sation of the unstable phonon mode<21 The displacements 
can be inferred from the corresponding eigenvectors and 
are sketched in Fig. [3 They belong to motions of the Fe- 
atoms either in direction of one of the two nearest neigh- 
bor Ni pairs, or perpendicular to this direction. If one 
describes the primitive cell as an Fe-octahedron encaged 
by Ni-atoms at the corner positions of the cube, the mode 
can be described in the first case as a rotational mo- 



tion of the perimeter atoms in the (001) Fe-plane, while 
the top and bottom atoms in the mixed plane remain in 
place. This kind of rotational mode is well known from 
perovskite structures, where oxygen octahedra form such 
tilted structures . 98 i" In the second case, it is rather de- 
scribed as a kind of breathing motion of the perimeter 
atoms in the Fe-plane, where one half of the atoms move 
inwards and the other half move outwards leading to an 
orthorhombic distortion of the Fe-octahedron. Keeping 
the terminology of Noda and Endoh as well as Kastner 
and coworkers, we will call these modes M4 and M2, 
respectively. Displacements according to both phonon 
modes lower the symmetry of the crystal and reduce 
the space group of the cubic LI2 phase, Pm3m (num- 
ber 221), to P4/mbm (number 127). Since the latter is 
already tetragonal, a variation of c/a does not imply an- 
other change of symmetry, which corresponds well to the 
experimental reports of a low temperature fct phase. 

As mentioned above, the nature of the distortions lead- 
ing to the softening of the phonon modes can be in- 
ferred in principle from the eigenvectors of the dynam- 
ical matrix which can be constructed, e.g., by a Born- 
von-Karmann fit to the experimental dispersions. This 
indirect procedure often leaves room for ambiguities and 
has lead to a controversy about the nature of the un- 
stable mode in the pasti 24 ' 27 A more direct approach is 
to monitor the variation of the total energy as a func- 
tion of the displacement according to the (frozen) phonon 
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modes. This can be realized in a straight forward manner 
within a density functional theory approach. We con- 
struct 2x2x1 supcrcells by doubling the primitive LI2 
cell in the plane of the atomic motion and apply the cor- 
responding displacements S as sketched in Fig. [7] The 
result, shown in Fig. unambigously proves that for all 
three alloys the M2 mode leads to a lower energy even 
for very small <5, explaining the imaginary frequencies at 
the M-point. In FesPt and Fe 3 Pd, the energy gain can 
reach values of about 5meV/atom for displacements of 
about 0.08 to 0.09 A. In Fe 3 Ni this energy gam is sig- 
nificantly smaller and the minimum energy is obtained 
for smaller displacements. On the other hand, the M4 
mode is considerably softer in FesNi than in the other 
two alloys. Here, the two lowest frequencies correspond 
to the two pure Fe modes, while for the other two al- 
loys the second eigenvalue has a marked contribution of 
the heavier element. In the latter case, the M4 mode is 
represented by the 5th and the 6th eigenvalue with fre- 
quencies of 18.7 meV and 19.4 meV for FesPd and Fe 3 Pt, 
respectively, as opposed to 9.3 meV for FesNi. 

In this respect, it may be important to notice that 
a translation of the mixed (001) plane in [110] direction, 
which leads to an exchange of the Fe and Ni-group atoms, 
transforms both modes, M2 and M4, in each other. Ap- 
plying this transformation once every second unit cell 
transforms the Ll 2 in DO22 order, as can be seen from 
Fig.[TJ DO22 can thus also be regarded as an LI2 structure 
with maximum density of anti-phase boundaries. Thus 
atoms in the Fc-plane will feel a superposition between 
both environtments and displacements according to M2 
and M4 modes become equivalent. Therefore, we expect 
that such distortions of the order will likely harden the 
soft phonon at the M-point, especially for FeaPd and 
Fe 3 Pt. 



C. Electronic origin of the antiferro-distortive 
transformation 

While the total energy calculations supply us with the 
information that the imaginary phonon mode is related 
to an orthorhombic distortion of the Fe octahedra in the 
LI2 unit cell, we are still left with the task to find the 
origin of this instability. A clue can be obtained from the 
comparison of the respective electronic density of states, 
which is shown paradigmatically in Fig. [9] for the undis- 
torted, perfect Ll 2 structure of FeaPt as well as for super- 
cells with M2 and M4 distortions. For all three cases, the 
majority spin density of states remains similar. The LI2 
DOS is in excellent agreement with the previous exten- 
sive study of Podgorny^ apart from a small shift of the 
Fermi level in both spin channels which can be attributed 
to the use of an exchange correlation functional without 
gradient corrections in Ref. [73|. FeaPt is at the onset of 
strong fcrromagnctism, thus the d-statcs of the majority 
spin channel arc filled and the respective DOS is small at 
the Fermi energy, Ep, leaving no room for changes that 



could account for the structural distortions. This, how- 
ever, does not apply for the minority spin channel, where 
a considerable number of e?-states is encountered right 
at E-p, especially for the undistorted LI2 structure. The 
vast majority of these states comes from the Fe atoms, 
while the Pt DOS remains small, due to the nearly filled 
rf-shcll and only moderate induced magnetic polarization 
of 0.36 ^b- 

On the other hand, the condensation of the M2-modc 
opens a deep pseudo-gap right at the Fermi-level which 
is the indication of an extensive redistribution of states. 
In the spirit of a band- Jahn- Teller mechanism, lifting the 
degeneracy of electronic states at the Fermi level can lead 
to a net gain in band energy. This has been previously 
proposed as the main electronic mechanism leading to 
the tetragonal distortion in disordered magnetic shape 
memory Fe-Pd alloys^ 

Accordingly, a distortion corresponding to the M4 
mode does not decrease the number of states at Ep. Nev- 
ertheless, it opens another pseudo-gap 0.25 eV below the 
Fermi level. Although this gap is very narrow and should 
be expected to vanish already at rather low temperatures, 
it appears possible that a moderate decrease of the va- 
lence electron concentration e/a may bring this pseudo- 
gap to the Fermi level and possibly causes a change the 
order of the modes at the M-point which alters the cor- 
responding ground state structure. 

The appearance of a pseudo-gap at Ep implies a se- 
vere reconstruction of the Fermi surface. Such recon- 
structions in connection with anomalous softening of the 
acoustic modes arc a sign of strong electron-phonon cou- 
pling and frequently related to a Kohn anomaly. This 
involves a significant number of nesting states at Ep that 
can be connected with the same reciprocal vector which 
describes the anomaly in the phonon dispersions. Fermi- 
surfacc nesting and strong electron-phonon coupling in 
combination with transitions to modulated phases are 
commonly found in conventional and magnetic shape 
memory systems, as Ni-Al, Ni-Ti and Ni-Mn-Ga^' 100 ' 101 
In fact, the phonon softening in disordered Fe-Ni Invar 
has been related previously to a Kohn anomaly in the 
nonmagnetic state on the basis of first principles calcula- 
tions of ordered Fe3Ni4S Here, we consider spin-polarized 
minority spin Fermi surfaces, i.e., the isosurfaces formed 
by the Kohn-Sham eigenvalues at the Fermi level in re- 
ciprocal space. These are depicted for FeaNi and FcaPt 
in Fig. 1101 and we find good agreement of the shape with 
previous results obtained using Wien2Kf^22 

Six bands are crossing Ep in the minority spin channel 
making the Fermi surfaces quite complex objects. Nev- 
ertheless, it is possible to relate nesting behavior to two 
of these bands, the 13th and the 15th band which are 
shown separately in Fig. 111! The 13th minority spin band 
represents an extended spherical hole pocket around the 
T-point with hopper-shaped extrusions reaching the Bril- 
loin zone boundary. Facilitated by the cubic symmetry, 
the flattened, horizontal walls of these extrusions can pro- 
vide nesting with a vector of exactly (i, ^,0), c.f., Fig. 




FIG. 9: (color online) Comparison of the total and partial, site resolved electronic density of states (DOS) of FesPt with perfect 
LI2 order (left) and M2 distorted FesPt with an displacement amplitude 8 — 0.795 A (center) and M4 distorted FeaPt with 
S = 0.748 A (right). The majority spin DOS are denoted by positive values, while negative values refer to the minority spin 
channel. In the M2 case, a pseudo-gap opens at the Fermi level due to a redistribution of the Fe states. For the M4 distortion, a 
narrower gap opens up 0.25 eV below the Fermi level, while the DOS at the Fermi level remains largely unchanged (pseudo-gaps 
denoted by arrows). 



[TT1 (a). These parts are connected by a tangential nest- 
ing vector, which is insensitive to a specific value of q 
and not expected to produce a sharp contribution to the 
generalized susceptibility. There is another possibility 
for nesting between the Fermi surfaces of 13th and the 
15th band. The 15th band forms a cubic electron pocket 
around the R-point as depicted in Fig. [TT](b). The walls 
of the cubes are nearly located on the same plane as the 
flattened parts of the 13th band. This plane is parallel to 
the boundary of the cubic Brillouin zone. Two of them 
exist in each of the cubic directions, such that nesting 
can occur between the two bands on the two different 
planes, respectively. The correspondence in this case is 
not necessarily exactly (i,i,0) and in direction of the 
face diagonal, but it can be expected to provide a pro- 
nounced peak in the generalized susceptibility and thus 
a considerable contribution to electron-phonon coupling. 
The cross sections in Fig. [TT] (c), taken at two differ- 
ent values of k z , demonstrate that this correspondence 
is not restricted to the zone boundary, but rather affects 
a significant portion of the Fermi suface. For FesNi the 
hopper-shaped extrusions of the 13th band are smaller 
and the cubes larger than in the Fermi-surface of FeaPt, 
which results in a smaller distance between the parallel 
planes and thus a smaller length of the q vector. This 
coincides nicely with the shift of the instability in the 
TAi -branch away from the M-point as well as the smaller 
magnitude of the imaginary eigenvalue. 



It should be mentioned that a thorough investigation 
of the relation between Fermi surface nesting and trans- 
formative processes requires a detailed statistics, e. g., an 
evaluation of the generalized susceptibility with respect 
to all contributing reciprocal vectors as well as the calcu- 
lation of the electron-phonon coupling matrix elements. 
This is a formidable task in itself and beyond the scope of 
this work. Nevertheless, in our view, the purely graphical 
discussion above gives already sufficient arguments that 
a Kohn anomaly is responsible for the low temperature 
structural changes in the ordered alloys. 



V. CONCLUSIONS 

Within the present investigation, we systematically 
compared structural, electronic and dynamic properties 
of ordered Fe-rich alloys with the elements of the Ni- 
group. The only one of these three alloy, which can be 
obtained experimentally in the ordered state, FesPt, is 
also the only one which possesses an fee ground-state 
structure. The dominant type of ordering in all three 
alloys is LI2, in accordance with experimental observa- 
tion. Competing DO22 order was observed for hypotheti- 
cal FeaNi and FeaPd for fct structures with a tetragonal 
distortion c/a > 1 and for FeaPt with bec coordination. 
However, for ideal fee (c/a = 1), DO22 order is well within 
the range of thermal energies, suggesting that anti-phase 
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FIG. 10: Plots of the isoenergy-surface in reciprocal space of the minority spin electronic states the Fermi level (Fermi surface) 
of FeaNi (left) and FeaPt (right). The Fermi level is intersected by 6 bands (shown online in different colors). The features are 
largely similar and differ mainly only in size due to the varying atomic volumes. Both surfaces exhibit large nearly flat portions 
which may give rise to a Kohn anomaly. 



boundaries are very likely to appear in the cubic LI2 
phases. 

The phonon dispersions of the three alloys possess as 
a common feature a complete softening of TAi branch in 
[110] direction around the M-point, which is strongest for 
FesPt and FesPd as well as a partial softening along the 
[111] direction. The first leads to the freezing of the re- 
spective phonon mode (M2) which can be identified as an 
orthorhombic distortion of the Fe octahedra in the prim- 
itive LI2 cell without involving the Ni-group elements. 
The second anomaly becomes increasingly pronounced 
with decreasing mass of the Ni-group element, as this 
leads effectively to an increased Fe-contribution to the 
low energy acoustical phonons. 

It may be speculated whether the rather smeared out 
softening of the TAi branch in [110] direction in the dis- 
ordered alloys can be interpreted in terms of a corre- 
sponding distortion with similar electronic origin, but on 
a larger, statistically modified period. This owes to the 
fact that the Fe-octahedra can be considered to be sta- 
tistically distributed within the bulk alloy and result fi- 
nally in a rather broadened and less pronounced anomaly. 
These octahedra are characteristic for the LI2 structure, 
but also other locally Fe-enriched cluster configurations, 
which arc more or less susceptible to deformations, may 
contribute. Similarly, the gross electronic structure will 
appear rather smeared out for the disordered crystal, but 



distinct features may prevail locally within suitable Fe- 
clusters and induce corresponding local distortions. In- 
deed, characteristic relaxations of the atomic positions 
have been reported in ab initio calculations for fee FeNi- 
Invar alloys for 64-atom quasi-random structures This 
also fits well with our observation that in FesNi not only 
the M2 mode but also the orthogonal M4 mode becomes 
rather soft at the M-point, which makes collective mo- 
tions in the Fe-(001) planes very cost effective. Fur- 
thermore, own investigations considering structural re- 
laxations of 108-atom Fe-Pd and Fe-Pt supercells with 
random distribution of atoms, yield a considerable en- 
ergy gain from a similar relaxation process^! 

The sum of these observations supports the view that 
the origin of the anomalous features in the acoustical 
branches can be attributed solely to the Fe species. This 
is confirmed by the calculation of the electronic struc- 
ture of the LI 2 ordered and the distorted alloys. The 
condensation of the M2 mode opens a pseudo-gap in the 
minority density of states right at the Fermi level, re- 
ducing the density of states by nearly 50%. Again, only 
contributions from Fe states are involved. The origin for 
this strong reconstruction of the Fermi surface can be 
related to extended flat areas which nest with a recip- 
rocal vector of nearly (^,4,0) which corresponds to the 
softening of the TAi phonon mode at the M-point. 

The subtle interdependence between electronic struc- 
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FIG. 11: Fermi surfaces of the 13th (a) and 15th (b) minority spin band of Fe3Pt (with an additional semitransparent image 
of the 13th band in the latter case). These two bands exhibit flat parts which may be connected with a wave vector close 
to (5,^,0) as indicated by the arrows. The vector can either be applied tangentially between the flattened, horizontal walls 
of the cone-shaped extensions at the zone boundary entirely within the surface of the 13th band as shown in (a) or between 
the vertical walls on the left and the respective sections of the cube shaped surface centered around the R-point in (b). The 
eigenvalues in reciprocal space of the of the 13th and 15th bands along two horizontal planes (k z =0.5 and k z =0.43) are shown 
in (c) as contour plots. The intersection of the complete Fermi surface (all bands) with these planes are marked by black lines. 
The arrows denote again parallel vectors (^,|,0). The energies are given relative to the Fermi level; the difference between two 
contour lines is 50 meV. 



ture at the Fermi level and the antiferro-distortive trans- 
formation confirms once more that the valence electron 
density, which controls the overall filling of the minority 
spin channel on a fine level, is one of the key factors for 
the design of Fe-rich functional alloys. On the basis of 
our investigation, we further conclude - in accordance 
with Ref. [7a - that the specific concentration of the Fe- 
species must be considered as another important quan- 
tity and, additionally, the mass of the alloying elements 
which modifies the lattice dynamics related to the low 
lying acoustical modes. 
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